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ABSTRACT 

It is well known that the combination of barotropic rotation and radiative equilibrium are mu¬ 
tually incompatible in stars. The Sun’s internal rotation is far from barotropic, however, which 
allows at least the theoretical possibility that the Sun’s thermal balance is one of radiative equi¬ 
librium in the region of the tachocline near the outer boundary of the radiative zone. We show 
here that (i) the constraint of radiative equilibrium leads to a straightforward ordinary differ¬ 
ential equation for the Sun’s rotation prohle, and (ii) solutions of this equation can be found 
that, to within current levels of accuracy, closely resemble the rotation prohle deduced from 
helioseismology. More generally, we calculate how large a baroclinic deviation from uniform 
rotation is required to maintain radiative equilibrium without meridional circulation through¬ 
out the bulk of the radiative zone. Very little deviation is required, well below detectability. 
The feasibility of radiative equilibrium for the tachocline suggests that the issue of a spreading 
boundary layer may be less severe than previously thought. 

Key words: hydrodynamics - radiative transfer - Sun: helioseismology - Sun: interior - Sun: 
rotation 


1 INTRODUCTION 

Since the availability of helioseismological data has allowed the re¬ 
construction of the details of the internal rotation of the Sun, signif¬ 
icant effort has been directed towards understanding what sustains 
the observed rotation rate Q.(r,9). The shape of the iso-rotation 
curves in the bulk of the convective zone was found to be well de¬ 
scribed by the characteristic solutions of the thermal wind equation 
( [Balbus et al.|200^ . In contrast, the physics of the relatively thin 
transition between the radiative and convective zone, known as the 
solar tachocline, is still largely uncertain. 

The well known theorem of Von Zeipel states that a star 
in uniform rotation cannot be in radiative equilibrium (see e.g. 
|Schwarzschild| ( |T958[ l). This is generally applied well within radia¬ 
tive cores, where any resdidual thermal energy imbalance is thought 
to be compensated by means of a mean velocity flow known as 
Eddington-Sweet circulation. However, [Balbus & Schaan] ( [2012[ ( 
have noted that matters are not so simple near the outer boundary of 
the radiative zone in a Sun-like star. The steady entropy equation, 
allowing for advection and radiative diffusion, is: 

P(u-V)cr =-(7-1)V-F, (1) 

where P is the pressure, u is the circulation velocity, a = 
log(Pp“^) is the entropy variable, 7 is the adiabatic index and 
F is the radiative flux. The radial, dominating component of Vcr 
is positive in the radiative zone and negative in the convective zone 
as for the Schwarzschild criterion, so that it has to go through zero 
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at the radiative-convective boundary. Therefore, there are two pos¬ 
sibilities at the boundary surface: either the circulation velocity is 
high enough for the (u ■ V)(7 term to be relevant and the angu¬ 
lar component of Vcr is finite, or the (u ■ V)(t term is negliglible 
and equation Q reduces to V ■ P = 0, the condition of radiative 
equilibrium. Neither of these possibilities is consistent with uniform 
rotation (or barotropic rotation more generally), and indeed the Sun 
shows a strong departure from rotation on cylinders in this region 
(and elsewhere). 

The Eddington-Sweet circulation velocity in the radiative zone 
of stars has a prominent role in the literature for solving the problem 
of energy transfer under conditions of uniform rotation. However, 
the inclusion of meridional circulation in the equations of stellar 
structure generates further complications (see [Tassoul[ ( [2000) for a 
discussion). The problem originates from the time-steady azimuthal 
component of the Euler equation of motion: 

(v ■ V)v =—IvP + g, (2) 

P 

that, for V = QR<j) + u, where H is the angular velocity, R the 
cylindrical radius and <j) the azimuthal direction, gives 

u-V(HP®) = 0, (3) 

i.e. for a star in uniform rotation, the circulation is along cylindrical 
surfaces at fixed distance from the rotation axis. This is incompat¬ 
ible with a model that is steady and symmetric with respect to the 
equatorial plane, as material would accumulate onto such a plane 
due to the circulation velocity. One is then compelled to consider 
non-steady models or to introduce an accommodating magnetic 
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field structure that would change equation Building a simple, 
comprehensive model for rotating stars when meridional circula¬ 
tion is included is not straightforward. 

For these reasons, it is interesting to ask whether the current 
data on the rotation of the outer part of the radiative zone of the Sun 
are compatible with a model in strict radiative equilibrium. Even in 
the bulk of the radiative zone, which is in near uniform rotation, the 
current accuracy of the data certainly does not allow one to say that 
it is rotating exactly as a rigid body. A small amount of differen¬ 
tial rotation could be present. Would this undetected deviation from 
uniformity be enough for the Sun to be in radiative equilibrium? 

The accuracy of the data of helioseismology decreases with 
depth in the Sun, and it is particularly difficult to study the prop¬ 
erties of a narrow transition layer like the tachocline, whose thick¬ 
ness may appear to be broader than it actually is because of the 
resolution of the inversion techniques involved in the data analysis. 
There are several different estimates of the width of the tachocline 
(see |Christensen-Dalsgaard & Thompson| p007[l), all i n the range 
1% —5% Rq; the most recent analysis by |Antia & Basu] ( |201 l| l sug¬ 
gests the possibility of a significant dependence of the tachocline 
width on solar latitude, with an average value not larger than 3% 
Rq but a strong uncertainty for the value at mid-latitudes (see their 
figure 2). This is thinner than the transition apparent from the data 
used in this paper. It is likely that the structure of the tachocline 
is still substantially unresolved. Any direct comparison of theory 
with the data must therefore be considered as a test of the plausibil¬ 
ity of the model, not as a precision attempt to reproduce the actual 
n(i?, z) in the interior of the Sun. The uncertainties are too signifi¬ 
cant. 

A discussion of the current knowledge of the internal rotation 
of the Sun through the radiative zone to the core is given in |Howe| 
( |2009| >. Estimating the systematic error in the rotation data for the 
interior of the Sun is not easy because of the complex inversion 
procedure involved jChaplin et al.|[T99^ . It is currently accepted 
that the data for the bulk of the convective zone of the Sun are very 
accurate, while those for the interior core are very uncertain. Al¬ 
though a thorough estimate of the errors in the outer radiative zone 
is not available, it is likely to be intermediate between the convec¬ 
tive zone and the core, with a relative uncertainty that is perhaps of 
order 10%; if one is conservative, possibly more (R. Howe, private 
communication). 

With these caveats in mind, in this paper, we construct a model 
of the region of the Sun interior to the convective zone under the as¬ 
sumption of radiative equilibrium, use this constraint to compute 
the rotation profile itself, and compare the resulting angular veloc¬ 
ity curves with the observed differential rotation of the Sun. Signif¬ 
icantly, we find that this model is compatible with the current data. 


1.1 The spreading of the tachocline 

An influential analysis of the tachocline was conducted by|Spieg^ 
|& Zahn| ^1992^ (see also |Zahn|2007t at a time when the early re¬ 
sults of helioseismology showed that the rotation of the Sun turned 
from strongly differential to uniform in a thin, unresolved transition 
at the radiative - convective boundary. These authors considered 
the standard equations of stellar structure, including the effect of 
circulation velocity. A time-dependent, initial-value problem with 
boundary conditions at the tachocline emerged. In the Spiegel & 
Zahn calculation, the term that offset the advection of angular mo¬ 
mentum in equation 0 above was dtl/dt, the time-explicit inertial 
derivative in the azimuthal equation of motion. Time-dependence 
thus played an essential role in the analysis from the start. 


The inclusion of this term led to the conclusion that the 
tachocline transition must spread with time. Specifically, [Spieg^ 
|& Zahn| \1992) showed that under these assumptions, a time- 
dependent diffusion-like equation for the angular velocity followed, 
implying ongoing penetration of the differential rotation into the ra¬ 
diative core of the Sun on a time-scale of r ~ 10® yr, significantly 
shorter than the age of the Sun. This posed a fundamental problem: 
why is the tachocline so thin? In the last few years, a significant 
amount of theoretical work has been devoted to studying candi¬ 
date mechanisms that would confine the spreading of an infinites¬ 
imally thin tachocline, both purely hydrodynamical (reviewed in 
Hughes et al.|200'7|l and magnetohydrodynamical jGough & Mcln-| 


tyre|1998 Garaud|2007^. However, the helioseismology data have 


not improved enough either to resolve the tachocline or determine 
a lower boundary for its thickness, or to significantly constrain its 
properties. 

The differential rotation of the Sun has been examined by 
numerical hydrodynamic and magnetohydrodynamic solar models, 
and important results have been obtained in reproducing rotation 
profiles for the convective zone that qualitatively resemble the re¬ 
sults of helioseismology. |Miesch et al.| j2006[ l imposed ad-hoc en¬ 
tropy boundary conditions at the base of the convective zone, and 
more recently there have been 2D (|Rogers|201 If and 3D jBrun et al.| 
|2011^ simulations of the coupled convective and radiative regions. 
The simulations are limited by the necessity of including artificially 
large diffusivities, which prevent an accurate force balance in the in¬ 
terior and in the tachocline, and coping with a wide range of time- 
scales, from those of convection and internal gravity waves to the 
predicted long time-scale of tachocline spreading. The results of 
|Rogers|(2011[ | (who finds that hydrodynamic processes can signif¬ 
icantly slow down the spreading of the tachocline) and |Brun et al.| 
POI 1[ > (who find a very fast spreading of the tachocline but relate it 
to a large viscous diffusivity), while significant on their own terms, 
may not fully capture the true dynamics of the radiative-convective 
boundary. 

The result by |Spiegel & Zahn] ( |1992| l implies that a thin (nearly 
discontinuous) tachocline, with initial and boundary conditions for 
the base of the convective zone drawn from helioseismology, is not 
a steady solution for the Sun, but instead gives rise to “burrow¬ 
ing” by differential rotation. This does not rule out the existence of 
a model of a steady, thin tachocline with a self-consistent angular 
velocity differing from the |Spiegel & Zahn|jl992^ profile, but nev¬ 
ertheless compatible with the present helioseismology data. In this 
paper, we present such a model. 

An outline of the paper is as follows. In §2 we present a de¬ 
tailed analysis of the governing equations, arriving at a set of two 
coupled ordinary differential equations that, with suitable boundary 
conditions, allow the angular velocity to be reconstructed. In §3, 
we apply our results to the solar tachocline and radiative interior. 
Einally, in §4 we summarize our conclusions. 


2 DESCRIPTION OF THE MODEL 
2.1 Equations of stellar equilibrium 

The steady equations of stellar equilibrium in the outer radiative 
zone for an inviscid flow are given by: 


(t) ■ Vlv = — - VP -b g, 

P 

(4) 

P- ^ keT, 

(5) 


fimp 
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V • Frad = 0, 

(6) 

Frad = VT, 

6pK 

(7) 


where v is the rotation velocity, p the density, P the gas pressure, g 
the (self) gravitational field, the atomic mass unit, p the mean 
molecular mass, T the temperature, Frad the radiative flux, a the 
radiation constant, c the speed of light, and k the Rosseland mean 
opacity. For a given chemical composition k is a function of the 
local density and temperature: k = k{p, T). 

We denote the rotational velocity in spherical coordinates 
(r, 6,(f))hy V = Qr sin 9(j), introducing the angular velocity D. The 
rotation is considered to be a small perturbation to a static spheri¬ 
cal structure, so that the equations can be linearized in the Eulerian 
perturbation variables 5P, Sp, Sg. 

sin^ 9 — - — -Sp + Sgr = 0, (8) 

p or p 

fl^rsin^cosS -— + 8gg = 0, (9) 

rp 69 


if the perturbations of the structural variables &P and have terms 
up to n -f 2 in their expansion: 


SP(r,9) = 

Pn{r) cos" 9, 

n=0,2,4 

(14) 

Sp{r, 9) = 

p„{r) cos" 9 . 

(15) 


n=0,2,4 


2.3 Solutions of the model in radiative equilibrium 

Expressions (T^l, i[T4]i and GU can be substituted into equations 
to (in which iT is expressed in terms of 5P, Sp by means 
of ITTJ) and the sum of all the terms proportional to the same power 
of cos 9 must be zero. Equation then yields for n > 2: 

P-n{r) = -ir^pn^_ 2 (r), (16) 

and equation ID gives, again for n > 2: 

Pn(r) = ^ ^pr(fl^(r) - n^_ 2 (r)) - • (17) 




ST SP Sp 
^ ^ "P~^’ 


( 11 ) 


where the first two equations are the r and 6 components of the Eu¬ 
ler force balance equation. The unsubscripted p, P and T quantities 
refer now to the unperturbed, nonrotating spherical equilibrium so¬ 
lution, and the perturbed Rosseland opacity Sk is given by: 

In the following analysis, we will neglect the self-gravity pertur¬ 
bations Spr, Sgg in equations ID and following the Cowling 
approximation. The validity of thie approximation is discussed in 
in Appendix A. 


2.2 Expansion in cos 6 

To solve equations ID-IH) in the variables Sp, SP and ST, it is 
convenient to expand these quantities in even powers of cos 9. The 
solar angular velocity given by the GONG data ( |Flill & et al.|1996] ) 
in the outer part of the radiative zone is very well approximated by 
a function of the form: 

0^(r, P) = Oo(f) + fl2(r) cos^ 0. (13) 

The accuracy of this fit may be seen in figures[T]and|^ which show 
the angular velocity isocontours of the Sun according to the GONG 
data (left) and an expansion of the form GD for r < O.767?0 in 
which Oo(f)> (^2(7) are 5-th order polynomials in r optimised to 
fit the data (right). The functions tlo{r) and are shown as 

the solid lines in figures We have used an expansion of the 
type l |13^ for the model of the radiatively sustained angular velocity 
curves. 

It has been noted by |Balbus & Schaan| ( |2012| l (see their equa¬ 
tion (85)) that if (r, 9) can be expressed in terms of the first n 
powers of cos^ 9, the equations of equilibrium can only be satisfied 


The perturbations in the pressure and density can therefore be ex¬ 
pressed in terms of the functions Gn(r) from the expansion of 
td{r, 9), except for the spherically symmetric terms Pair), Po(j'). 
Knowledge of those terms is not required to impose radiative equi¬ 
librium, and they can ultimately be considered as absorbed into the 
non rotating, spherically symmetric solution. 

We may eliminate from the above, obtaining a relation be¬ 
tween pn and Pn : 

1 dP 

rpn = - [nP-n - (nP 2)P„+2]-(18) 

r dr 


The same technique applied to the radiative equilibrium equation 
GD leads to the rather cumbersome form 

+ 1 d ( r^T'^ V _ (Pn 


pr'^ dr 
C d 

J .2 

rp4 


Pr^ dr\P 


r^pK 


(7) 

(n-h 1) [(n-f 2) 


■2 dr \ 


Ok , _ 9k / Pn p- 
Pn^ +T - (- 


dT\ P 
Pn + 2 


P 

Pn+2 


pn 

+ 


m-m- 


)]} 


= 0, 

where the constant 


(19) 


SLq 

167rac ’ 


( 20 ) 


has been introduced, where Lq is the Solar luminosity, and we 
have used prime notation [X)' for dX/dr to aid readability. This 
equation holds for n = 2,4, with the understanding that Pe = 0, 
P6 = 0. Equations G§ and 1^, with proper boundary conditions, 
are coupled ordinary differential equations that uniquely determine 
the solution of our problem. In practice, we have found that the best 
way to proceed is to express Pn{r), pn{r) in terms of Go(f)> 1^2 (r) 
via equations GD and ([TtJ and use equation l |19^ for n = 2,4 to 
obtain a set of two lengthy, coupled ordinary differential equations 
for the This is a straightforward procedure, but we will not ex¬ 
plicitly write out the lengthy equations here, as there is little to be 
gained beyond the content of equations l |18[ l and l |19^ . The equa¬ 
tions may be solved, once boundary conditions are specified for Gq, 
GI and their first and second derivatives at some radius ro. 
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GONG DATA 





GONG DATA & FIT 
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Figure 1. Isocontours of the angular velocity f2(r, 9) for 0.55i?Q < r < Rq. The left panel shows the GONG data. The right panel shows our best fit of the 
form (ID for the data interior to r < 0.76i?Q, joined with the actual GONG data for r > 0.76Rq (the surface r = O.IQRq is shown as a black line). The 
scale is in rad s~^. Figure[^shows a more detailed comparison of the values in the interval 0.557^0 < r < O.76i?0. 



Figure 2. Isocontours of the angular velocity 0(r, 9) for 0.56Rq < r < O.76i?0. As in figure^ the left panel shows the GONG data, while the right panel 
shows our best fit of the form (ID for the data. The scale is in rad s ^. 


3 APPLICATION TO THE OUTER RADIATIVE ZONE 

In the following analysis, the equilibrium variables P, p and g have 
been taken from the solar model of |Bahcall et al.| ( [2005| l. The Rosse- 
land opacity k{p, T) and its derivatives have been interpolated from 
the OPAL table for solar composition Pglesias & Rogers|1996^ . 

Our approach is to solve the coupled tP equations for compar¬ 
ison with the rotational GONG data. We set the boundary conditions 
for Q? and its first and second radial derivatives at ro = 0.607?©, 
where the rotation rate is approximately uniform, and integrate out¬ 


wards. We adjust the spatial derivatives within the constraint of in¬ 
terior quasi-uniform rotation, and fit the curve to the GONG data. 
The fit has been executed by performing a minimization around 
an initial, exploratory solution. The resulting boundary conditions 
imposed at ro = 0.60 Rq are shown in table[T] The angular veloc¬ 
ity terms Oo(t), OKt) of equation l[^i for this radiative solution 
are shown as the dashed lines in figures along with the GONG 
data. 

Figures ^ indicate that solutions of the model in radiative 
equilibrium can be found which fit the data well and match closely 
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0.60 0.62 0.64 0.66 0.68 0.70 0.72 


Figure 3. Angular velocity term Qq (r) of equation (m in the outer part of the radiative zone according to the GONG data (solid line), to the model in radiative 
equilibrium that best fits the data (dashed line), and to models a (dots) and b (crosses); see description in the text. The boundary conditions at vq = 0.60 Rq 
ai‘e illustrated in table^ An en'or bar con'esponding to a relative uncertainty of about 5% at the tachocline is shown, although the actual en'or is larger. 



Figure 4. Angular velocity term (^) equation (T3 in the outer part of the radiative zone according to the GONG data (solid line), to the model in radiative 
equilibrium that best fits the data (dashed line), and to models a (dots) and b (crosses); see description in the text. The boundary conditions at vq = 0.60 Rq 
ai‘e illustrated in table^ An en'or bar con'esponding to a relative uncertainty of about 10% at the tachocline is shown, although the actual error is lai'ger. 
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9.x 10“^^; 
8.5x10“^^ 

8.x 10“^^ 

•r 7.5x10“^^ 

"c? 7.x10-^^ ’ 
6.5x10“^^ 

6.x 10“^^; 

5.5x10-^^ 

0.40 


- Radiative Baiance 


0.45 0.50 0.55 0.60 0.65 

r/Ro 

(a) 



(b) 


Figure 5. a\ Angular velocity $l^(r) in a Sun-like star, obtained imposing boundary conditions of r2Q(ro) = 7.5 • 10“^^ s~^ (a typical value of the angular 
velocity in the Sun), dO^^iro)/dr = 0, dP'O^^irQ)/dr'^ = 0 at ro = O.65i?0. b\ Deviation from uniform rotation expressed as (f7Q(r) — f2Q(ro))/(^^(ro)) 
with ro = O.65il0. 



Best fit 

Model a 

Model b 

Hg(ro) 

7.4 ■ 10“^^ 

7.4- 10“^^ 

7.4- 10“^^ 

dn‘f^(ro)/dr 

3.9 ■ 10-^^ 

2 . 0 - 10 -^^ 

7.8 ■ 10-^^ 

(i^A^(ro)/dr^ 

-4.5 ■ 10-“^ 

1 

0 

CO 

CM 

1 

-9.1 ■ 10“^^ 


-2.4- 10-1^ 

-2.4- 10“^^ 

-2.4 ■ 10“^^ 


-5.9- 10“^^ 

-3.0 ■ 10“^^ 

- 1.2 ■ 10 "^^ 

d?^^(rf) / dr^ 

- 1.2 ■ 10 “^^ 

- 0.6 ■ 10 “^^ 

-2.4 ■ 10"^^ 


Table 1. Boundary conditions imposed at ro = 0.60 R 0 for the solutions 
of figures[^|^ The units are in cgs and have been omitted to aid readability. 


to a uniformly rotating interior. We have shown illustrative 5% and 
10 % error bars for the values of ^2 in figures 1^1^ It is clear 
that the model fits well to the data within the expected uncertain¬ 
ties. This result raises the possibility that a very small deviation 
from uniform rotation might be sufficient to allow the deeper stellar 
interior {R < O.fiSi?©) to be in baroclinic radiative equilibrium. 

It should be noted that although the data from helioseismol¬ 
ogy impose a costraint on the values of flo(fo), they are 

too uncertain to constrain the values of the derivatives of fig. ^2 
deep in the radiative zone. It is therefore difficult to compare the 
derivative boundary conditions we impose with those derived from 
helioseismology. 

A test of the robustness and sensitivity of the solutions can 
be performed by noting how they change when the unconstrained 
boundary conditions drdj{ro)/dr, (PtlQ(ro)Jdr^, dtl^b'oi/dr, 
d^Q, 2 {ro)/dr^ are varied. We show in figuresl^Wthe solutions for 
the case in which the boundary conditions are a factor of 2 smaller 
(model a) and larger (model b) than those of the best fit. The respec¬ 
tive boundary conditions for each model are shown in table[T] From 
figures D and 0. it is clear that the validity of the assumption 
of radiative equilibrium is not acutely sensitive to possible future 
moderate revisions of the angular velocity profile. 

3.1 Deviation from uniform rotation 

As a preliminary investigation to determine what amount of differ¬ 
ential rotation is required for a star to be in radiative equilibrium, the 
model has been applied to a Solar-like star, i.e. a star with the same 
pressure, density and temperature profiles as the model by[Bahcallj 
|et al.| (|2005[> but a different rotation profile. We integrate towards 


f^o(^o) 

7.5 ■ 10 “^^ s-^ 

dA^{ro)/dr 

0 

d"ng(ro)/dr^ 

0 

Uf(rp) 

0 

/ dr 

0 


0 


Table 2. Boundary conditions imposed at rg = 0.65 R© for the solution of 
figurej^ 

the deep interior with boundary conditions at ro = O.fiS/?© locally 
resembling uniform rotation and using a value of fl(ro) close to the 
average angular rotation rate of the radiative zone of the Sun. The 
boundary conditions are shown in table|^ 

The model is able to maintain nearly uniform rotation down to 
a significant depth in the radiative zone by constraining (f) = 0 
(see figure]^. The solution starts to deviate more strongly (beyond 
the 1% level) only at a significant depth (with a deviation of 4% at 
r ~ 0.4 Rq). Were such a profile dynamically stable, the bulk of 
the radiative zone of a Sun-like star could be in radiative equilib¬ 
rium and quasi-uniform rotation, without the need for meridional 
circulation. The question of stability is further discussed below. 


4 CONCLUSION 

Meridional circulation is generally seen as a means to maintain 
thermal energy balance in rotating stars. However, the purpose of 
meridional circulation is to maintain uniform or barotropic rota¬ 
tion; if the rotation is baroclinic, it might well be compatible with 
radiative equilibrium. Nevertheless, small circulation velocities are 
usually invoked even in baroclinic flows to maintain thermal bal¬ 
ance. 

While solving one problem, meridional circulation seems to 
create another: starting with the influential study of |Spiegel & Zahn] 
\\992) , circulation-induced diffusive spreading of the tachocline 
has been viewed as a major problem for understanding solar rota¬ 
tion. In this paper, we have taken a step back by arguing that since 
baroclinic solar rotation is a reality, one should investigate what 
sort of rotation profiles emerge when the constraint of strict radia¬ 
tive equilibrium is applied, without the benefit (or complications) of 
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circulation velocities. We have shown that it is easy to find solutions 
that fit well to the observed rotation profiles. Indeed, it is not clear 
that even the deep radiative interior is free of baroclinic rotation at 
the small level needed to influence radiative balance. 

The most striking feature of the model presented in the current 
paper is its simplicity. It is devoid of meridional circulation currents, 
magnetic fields, viscosity and compositional gradients, effects of¬ 
ten invoked to explain the physics of the tachocline ( [Hughes et al.| 
[2007^ . It is by far the simplest interpretation that is compatible with 
the current data from helioseismology. Interestingly, however, the 
problem of time dependence might still implicitly be present in our 
models, since we have not addressed the all important question of 
the stability of our baroclinic rotation profiles. The stability of stel¬ 
lar rotation curves is complex because the criteria are significantly 
affected by such subtlelties as the presence of even weak magnetic 
fields jBalbus|1995^ , the thermal diffusion of the displaced fluid ele¬ 
ments I Goldreich & Schubert|1967), viscosity ( |Acheson|1978^ , and 
even resistivity jMenou et al.|20()^ . One advantage of our formula¬ 
tion is that it is simple enough to lend itself to a rigorous linear sta¬ 
bility analysis. If such an analysis shows that the profiles are stable, 
this would go some way to alleviating the problem of the spreading 
of the tachocline. A general three-dimensional, magnetised, nona- 
diabatic study of the stability of our solutions will be presented in a 
subsequent paper. 
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APPENDIX A: DISCUSSION OE SELF-GRAVITY 

In this appendix we estimate the 5gr, Sgg terms of equations 
and discuss the Cowling approximation. For this purpose, it will 
suffice to consider a star with uniform angular velocity If and to 
determine the perturbation in the gravitational potential (5<l> due to 
the rotation-induced structural change of the star. 

In this case it is more convenient to expand the pressure and 
density perturbations in terms of Legendre polynomials P; (cos 9) 
rather than in terms of powers of cos 9 as in equations CD,®: 

SP{r,9) ^ Po{r) + P 2 {r)V 2 {u), (Al) 


Sp{r, 9) = po{r) + p 2 (r)V 2 {p), (A2) 


with p = cos9. The Pair), P 2 {r), po(f), P 2 (f) terms of this ex¬ 
pansion differ from those of an expansion of the form jl4| (, jl5[ ( by 
a numerical factor. 

Plugging expressions | |A1| (, jA2[ ( into equations ([^, ne¬ 
glecting the Sgr, Sgg terms (hence the perturbation in the gravita¬ 
tional potential will not be self-consistent), and grouping the terms 
proportional to the same Legendre polynomials in each equation, 
it is possible to derive 7^2(r), P 2 {r). The procedure is analogous 
to that leading to equations ( [16^ , jl7[ l. The result for the density 
perturbation is: 


P2(r) 


1 dp 
3 g dr 


(A3) 


The perturbation to the gravitational potential in the star is: 

S^{r) = f (A4) 
J \r-r'\ 

Using equation l |A3^ and retaining only the non spherically sym¬ 
metric part, this is: 

iHr) = / y (AS) 

The integral can be performed by expanding: 

+ (A6) 


i>o > 


where the terms 

Q* = E A'), (A7) 

m^O 

do not contribute to the integral in l |A5| l. Plugging the expansion 
l |A6[ ( in equation ( |A5^ the integral is computed by making use of 
the Legendre polynomials orthogonality properties as: 


5 $ = - 


15 




r'® dp 
g dr' 


dr' + r^ 



r' dp 
g dr' 



(A8) 

By halting the integration at r' = Rq, we are neglecting the change 
of the shape of the surface of the star due to the rotation and the ef¬ 
fect of the material outside the surface of the non-rotating structure. 
This is justified as the density in the outer layers, and therefore their 
contribution to the perturbation in the potential, is low, as verified 
by our numerical computation of the integrals. 

We determined the potential perturbation (5<l> assuming the 
same Solar background structure as in sectionj^and a uniform an¬ 
gular velocity tA = 7.5 • 10“^^ s“^. The amplitudes of the result¬ 
ing gravity perturbations computed from the components of 
at r = 0.65 R© and p = 1 are: 


|(5c/r(O.657?0)| 


1^1 
I dr \ 


= 1.1-10 ^cm s 


(A9) 


I 1 dS^ I o o 

|5gr6i(O.65i?0)| = j -I = 1.5 • 10 ^cm s . (AlO) 

The ratios of these terms to the penultimate terms in equations l[^, 
are: 


I Sgr I 

I gSp/p\ 


= 1 %, 


(All) 


__I = 9% 

{l/{rp)){dSP/d9)\ 


(A12) 


We have therefore shown that while the radial self-gravity can 
safely be neglected, the horizontal one is small but not tiny. The 
Cowling approximation can be employed in the study of rotating 
stars, but precision work on their interior structure should be con¬ 
ducted retaining these additional terms. 
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